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We discuss an idea of how 3D critical exponents can be determined by Conformal 
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1 Introduction 

There are many reasons to study Conformal Field Theory (CFT): 

• Beyond the Standard Model physics can contain sectors which are conformal in the 
UV. These can be ad hoc sectors, like unparticles, or sectors having to do with the 
Electroweak Symmetry Breaking, like in Walking or Conformal Technicolor. 

• Understanding of quantum gravity in AdS can be achieved by asking questions about 
a dual CFT. Conformal symmetry also governs structure of field correlators in dS 
space, relevant for inflation. 

• Finally, conformal symmetry is crucial for the theory of critical phenomena, be that in 
classical statistical physics, or in quantum condensed matter ('quantum criticality'). 
This was the primordial reason for the birth of CFT. This is also the motivation behind 
this paper. 

Consider the paradigmatic model of criticality, the classical O(N) ferromagnet in D = 3 
dimensions, described by a partition function 



where o~i are unit iV-component vectors at the vertices of a regular cubic lattice. The 
simplest case N = 1, a = ±1, is the Ising model with the Z2 global symmetry. 

Near and slightly above the critical temperature, the spin-spin correlation function of 
this model is expected to be power-like at short distances with an exponential decay at large 
distances: 




(1.1) 




r >£, 



(1.2) 
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where £ = £(T) is the correlation length. At T — > T c the correlation length diverges, and 
the spin-spin correlator exhibits scale invariance. See [1] for this and other elementary facts 
about critical phenomena. 

The emergence of conformal symmetry at the critical point is more mysterious. This 
seems to be a generic feature of criticality but why this happens is not fully understood [2] . 
Recently there was a renewed interest in the question whether there exist interesting scale 
invariant but not conformally invariant systems [3-6]. At this time, we don't have anything 
to add to this discussion and will take conformal invariance of the critical point for granted. 



2 Critical exponents. . . 

The most important parameters characterizing the critical point are the critical exponents. 
In this talk we will mostly need only two of them. The first one is the spin field dimension A a 
defined in (1.2). 1 The second one is the energy field dimension A e , defined by an analogous 
equation 

( £ (r)e(0)) = r- 2A *. (2.1) 

The energy field e is interesting because it's the lowest-dimension scalar operator of the 
theory which is a singlet of the global symmetry (O(N) or Z2). This operator is relevant: 

A £ < D = 3 (2.2) 

(see Eq. (3.1) below). Adding it to the CFT action as a perturbation: 

Soft ->■ S CF T + ta A ^ 3 J d 3 xe(x), t<l, (2.3) 

we generate an RG flow. If t <C 1, then the perturbation is small at the UV cutoff length 
scale a, but it eventually becomes big in the IR, breaking conformal symmetry. The length 
scale £ at which this happens is determined from the equation 

£ A *- 3 ~ta A *- 3 . (2.4) 

This length coincides with the correlation length. 

In fact, Eq. (2.3) is the continuum space description of what happens in the Ising model 
when the temperature is slightly detuned away from the critical temperature, with t = 
(T — T c )/T c the reduced temperature. Eq. (2.4) gives then the well-known equation for the 
critical exponent u, which measures how fast £ grows when T — > T c , in terms of A e . 

It is also important that e is the only relevant singlet scalar. To reach the Ising critical 
point, we finetune just one parameter (the temperature), hence one relevant singlet scalar. 
More such scalars would imply multicritical behavior. 

lr This parameter is in one-to-one correspondence with the critical exponent commonly called r\. 
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3 ... and what is known about them 



All existing techniques to determine the 3D critical exponents are based on the idea of 
universality. Universality means that lots of models with different UV structure, once put 
at the critical point will reduce in the IR to the same CFT with the same critical exponents. 
The only thing which should matter is the global symmetry of the underlying Lagrangian 
(Z 2 , 0(N), etc.). 

One can take then any real-world system whose microscopic structure has the needed 
symmetry (magnets etc) and extract critical exponents experimentally, see e.g. [8] for a 
summary of existing measurements. 

In the same spirit, one can consider the lattice model (1.1) and do Monte-Carlo simula- 
tions. These "measurements" are computer time-consuming, but one has more control over 
interesting observables. One can also modify the lattice action to improve the numerical 
performance by reducing lattice artifacts etc. Using these techniques, Ref. [9] reports the 
following 3D Ising dimensions: 

A CT = 0.5183(4), A £ = 1.412(1) (3D Ising). (3.1) 

However, by far the most famous technique to compute the O(N) 3D critical exponents 
has been the e-expansion [10]. In this case the UV theory is the theory of N free scalars tp 1 
perturbed by the operator S = (</> 2 ) 2 , relevant if D < 4: 

^ = \{d^) 2 + Xa A ^S. (3.2) 

The IR theory in D = 3 should be the O(N) critical point, but to get there directly is difficult 
since the IR fixed point lies far away from the UV one in the theory space. Alternatively, 
one could say that the 3D O(N) critical point does not allow a weakly coupled description 
in terms of the UV degrees of freedom. 

The idea of Wilson and Fisher was that in 4 — e dimensions for e <^ 1 the O(N) critical 
point becomes weakly coupled. In this case the perturbing operator S is only weakly relevant: 

A s = (4-e)-e. (3.3) 

The lowest order /3-function takes the form 

/3 A = _ e A + (iV + 8)^, (3.4) 

and an IR fixed point lies at 

(3 ' 5) 

So, in 4 — e dimensions all quantities of interest (like the IR operator dimensions) can be 
computed perturbatively, expanding in e. For example, the first nontrivial terms for and 
A^2 are given by 

A^(l-e/2) + i |±| F e 2 + ..., 

V = (2-e) + ^±| C + ..., (3.6) 
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flow 




Figure 1: The accepted phase diagram of the O(N) fixed points for 3 < D < 4. 



where the terms in parentheses denote the UV dimensions. 

One then makes the natural assumption that this line of fixed points as a function of 
e continuously connects to the 3D O(N) fixed point in the limit e — > 1; see Fig. 1. In 
particular, the field should connect to a and 4> 2 to e. Setting e = 1 in the series (3.6) then 
allows to calculate 3D critical exponents. 

The series (3.6) have been extended to terms of order e 5 [8], and the 3D exponents 
evaluated this way are in a very good numerical agreement with experiment and Monte- 
Carlo studies. Notice however that as always with perturbation theory, these series are 
only asymptotic. For e = 1, the divergent nature of the series starts to show already after 
the first couple of terms, so they need to be resummed (via a Borel transform), otherwise 
higher order terms do not improve the accuracy. The resummed results have a residual 
theoretical ambiguity associated with free parameters entering the resummation procedure. 
To quote [8]: 

"Needless to say, with free parameters and short initial series it becomes possible to find 
occasionally some transformed series whose apparent convergence is deceptively good." 

4 Why one would like to do better 

The situation described in the previous section is not fully satisfactory for the following 
reason. Presumably, the O(N) critical points are very special isolated points in the space 
of theories. It would be nice to be able to pick out these "diamonds" by focussing on them 
directly, without having to flow from free theory. After all, this is what happens in the two 
dimensional case, where there is a neat classification of critical points based on conformal 
symmetry [1 1] . The critical exponents of the 2D Ising model can then be determined exactly: 



On the other hand, conformal symmetry is left unused in the RG calculations leading to 
the e-expansion. This is because this symmetry only emerges at the critical point; it's not 
present along the flow. 

I would like to quote from a 2003 interview with A.M. Polyakov [12], where the desire 
for a better 3D theory has been stated with clarity: 

"Let me tell you what I think of the renormalization group. I think there are two types of useful 
equations. One type is human-made, they are invented by people. The other type reflects some 
'pre-established harmony.' They can be discovered (uncovered) and not invented. Renormalization 



A CT = 1/8, A e = 1 




(4.1) 
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group is clearly a human made thing. It's clearly a smart way of calculating things but it doesn't 
have a breathtaking quality of, say, the Dirac equation. 

The example of the second kind is operator product expansions. They form some beautiful 
mathematical relations and I was dreaming in the 1970s to have some classification of fixed points 
based on the possible operator product expansions. The program was a little like classifying Lie 
algebras. In that case you start with the commutator relations which define the Lie algebra and 
then you classify all possible semi-simple algebras. You arrive at a stunningly beautiful theory 
(which was clearly discovered and not invented). I was working on that project in the 1970s and 
I still think it might have a chance. It was successful in two dimensions. We can classify possible 
fixed points in two dimensions using operator product expansions. That's what conformal field 
theories are about. And I think it's not excluded, that in 3 dimensions something like that is still 
possible. I was working for a while on this without much success in the 1970s and then I switched 
to other things. 

I think the epsilon expansion ended the subject in the practical sense. You can calculate more 
or less what you want with good accuracy but aesthetically the subject is not closed yet. It's 
possible that there will be classification of fixed points in three dimensions, based on string theory, 
similar to what we have in two dimensions. But that's just dreams." 

5 Conformal bootstrap 

We would now like to show how information about operator dimensions can be extracted by 
conformal field theory techniques. The idea is called "conformal bootstrap" and is actually 
quite old [13, 14]. However, its practical implementation has appeared only recently. In 
this section we will demonstrate results achieved so far in D = 4 and D = 2 dimensions. 
In the next section we will explain how the same idea applied in D = 3 could be used to 
determine the O(N) critical exponents. As we proceed, it should become clear why we hope 
to succeed where previous attempts did not have much luck. 

5.1 Old stuff 

Conformal bootstrap studies consistency conditions for four point functions. The simplest 
constraint comes from the four point function of the lowest dimension scalar, a in the Ising 
model case: 

(a(xi)a(x 2 )(r{x 3 )a(xi)) . (5.1) 
Conformal symmetry implies that this correlator must have the form 

x^ A °x^ A °g(u,v), (5.2) 

where x^- = \x~ Xj\ and g(u,v) is a function of the conformal cross ratios u = (x^ 2 Xg 4 ) / (^13^24), 

V = (^14^23) / '( x 13 x 24) ■ 

A representation for this correlator and for the function g(u, v) can be obtained by using 
the OPE. The leading terms in the a x o OPE are 

a(x)a(0) = \x\- 2A °(l + X aa£ \x\ A ^e(0) + ...), (5.3) 
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where the OPE coefficient X aa£ also appears as an overall constant in the three point function 

{a{x) a(0) e{y)) = R^^F^F ' (5 ' 4) 

whose coordinate dependence is fixed by the conformal symmetry [15]. 

Beyond the leading terms, the general conformally invariant OPE can be written as a 
sum over conformal primary operators & of the form: 

cr(xi) a(x 2 ) = ^ \fi C(xx - x 2 , d X2 )ff(x 2 ) . (5.5) 



In other words, each operator & appears in the OPE accompanied by its derivatives. The 
coefficient functions C{x\ — x 2 , d X2 ) giving these subleading contributions are fixed by con- 
formal symmetry. There is one overall coefficient A^ per conformal family. 

In general, the a x a OPE will contain infinitely many primaries = (?® of spins 
/ = 0, 2, 4 . . . Only even spins will appear since we are considering the OPE of two identical 
scalars. The dimensions of these operators and the coefficients A<? are at this stage free 
parameters. However, unitarity of the theory (= reflection positivity in the Euclidean) 
implies that \@ should be real, and that the dimensions of spin / fields should satisfy the 
lower bounds [16, 17]: 

A>{ D ' 2 - 1 <' = °>- (5.6) 
~ [l + D-2 (I > 1). 

Coming back to the four point function (5.1), one can apply OPE to the 12 and 34 pairs 
of points. One then reduces the four point functions to a sum of two point functions acted 
upon by the differential operators C{x\ — x 2 ,d X2 ) and C(x3 — x^d Xi ). The sum involves 
only diagonal terms: & — Iff' since only identical fields have nonzero two point functions. 
It follows that the function g(u,v) can be represented in the form 



g{u,v) 



^(X e ) 2 g^(u,v), (5.7) 



e 



where g^ are completely fixed functions of u, v which depend only on the dimension and 
spin of the exchanged operator iff. These functions are called conformal blocks. 

We are now ready to formulate the conformal bootstrap consistency condition: the func- 
tion g(u,v) represented as a sum of conformal blocks should satisfy the crossing symmetry 
relation: 



g(u,v) = (u/v) A °g(v,u) . (5.8) 



This relation follows from the fact that the four point function (5.1), (5. 2) should be invariant 
under 1 -H- 3, which corresponds to u •<-> v. 

Alternatively, the same condition can be formulated as the associativity of the OPE. We 
can do the OPE either in the 12-34 channel, or in the 13-24 channel, and the result should 
come out the same. 

A word of caution concerning terminology: It was probably Ref. [11] that first used the 
term "conformal bootstrap" in the sense just described (although the idea is 10 years older). 
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In more ancient times the same name was sometimes used to denote a completely different 
thing: Feynman diagram perturbation theory with bare propagators having anomalous 
dimensions built in from the start. These dimensions then had to be determined from a kind 
of Dyson- Schwinger equation. This technique is also known as the "skeleton expansion". 
It makes sense in the large iV approximation, see e.g. [18], but is not useful in general. To 
quote classics [14]: 

"The form of the [skeleton expansion] equations depended in an essential way on the type of 
fundamental fields and on the form of their bare interactions, whereas the results of the theory 
with anomalous dimensions should not be sensitive to the choice of the initial Hamiltonian" 

5.2 Recent stuff 

The bootstrap condition was first formulated in the 70s [13, 14] but at the time it did not 
lead to solutions of 3D or 4D CFTs. 2 It seems that two basic difficulties were responsible 
for the lack of progress. First, conformal blocks were not well understood at the time; they 
were known only in terms of complicated integrals or power series expansions in u, v. The 
second reason is that the crossing symmetry constraint (5.8) is objectively hard, being a 
functional equation for infinitely many operator dimensions and the OPE coefficients. It 
takes some effort to imagine how such an equation can be useful without further input. 
One may also fear that to get control over the theory one will have to also impose crossing 
for correlators of all fields rather than just cr, at which point the problem gets even more 
unwieldy. 

However, recently progress has been achieved both in finding simple expressions for 
conformal blocks and in extracting concrete results out of the bootstrap condition. 

First, Dolan and Osborn [19, 20] put the D = 4 conformal blocks in the following explicit 
form: 

9$M = ( -if—- [ W*)*A-i- a ® z)] (D = 4) , 

Zi % % 

u = zz, v = (l-z)(l-z). (5.9) 

Second, in [21] we studied constraints imposed by the bootstrap equation on A E as a 
function of A CT . That study was also performed in D = 4, motivated by some questions in 
electroweak phenomenology beyond the Standard Model, but this is not important for the 
present discussion. Rather surprisingly, we found a completely universal upper bound of 
the form 

A e </(A„), (5.10) 

starting at the 4D free scalar theory point (A CT = 1, A e = 2) and then growing monoton- 
ically; see Fig. 2. The bound was found numerically, using the explicit conformal blocks 

2 On the other hand, in 2D the bootstrap program was effectively carried out [If]. The simplification 
in 2D was that fiefds are grouped into bigger multiplets under the Virasoro algebra extending the finite- 
dimcnsionai conformal group. In minimal models, there is only a finite number of Virasoro primaries, hence 
a finite number of Virasoro conformal blocks. In addition the dimensions of the fields are known. For this 
reason crossing symmetry constraint could be solved, determining also the OPE coefficients. 
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free scalar 



D=4 



allowed region 



Figure 2: A bound on A £ as a function of A CT , shown here schematically. 



(5.9). Reality of \& (hence positivity of expansion coefficients in the conformal block rep- 
resentation) and the unitarity bounds (5.6) also played an important role in the analysis. 

Since then these bounds have been strengthened by improving the numerical techniques 
[22], extended to the case when a CFT has a global symmetry [26, 27], or supersymmetry 
[24, 27]. Also, the same techniques have been used to derive universal upper bounds on the 
OPE coefficients and lower bounds on the CFT central charges [23-25]. Further dramatic 
improvement of the algorithm has been recently achieved in [28], where current best possible 
bounds of all kinds can be found. 

The just mentioned bounds represent valid constraints on the landscape of possible 4D 
CFTs. However, it is not known at present if there are any nontrivial 4D CFTs saturat- 
ing them. The situation is different in 2D, where the corresponding bound has also been 
computed in [21], and then with better accuracy in [22]. This was possible because the 2D 
conformal blocks are even simpler than in 4D [19, 20] : 3 

g$(u,v) cc k A+l (z)k A ^(z) + (z^z) (D = 2) . (5.11) 



free scalar 



D=2 



Ising 



other minimal 
models 



4A CT 



Figure 3: Bound analogous to Fig. 2, but for D = 2. 



This 2D bound is given schematically in Fig. 3 and shows interesting structure. Just as 
in D = 4, the bounds begins at the free scalar point and grows monotonically. This growth 

3 These are of course the "small" SL(2, C) conformal blocks, appropriate for comparison with 4D. "Big" 
conformal blocks for full Virasoro algebra depend on the central charge and are much more complicated 
objects. 
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starts off linear rather than as a square root, tangent to the line 

A £ = 4A CT . (5.12) 

This fact has the following explanation [21, 22]. The free scalar theory in 2D contains 
primary operators V a = e iaX of dimension a 2 (for appropriately normalized X) and with a 
leading OPE 

V a xV a ~V 2a , V a xV£~l. (5.13) 
This implies an Ising-type OPE for the real parts of these operators: 

crxcr~l + e + ..., <j = W/ a , e = $tV 2a . (5.14) 

Since the so defined fields a and e have dimensions related by (5.12), this line must lie in 
the region allowed by our bound. That it is actually tangent to the bound means that our 
result becomes best possible in the limit A CT — > 0. 

However, the most interesting feature of the 2D bound is that a number of solvable 2D 
CFTs are found to saturate it. The first of these is the 2D Ising model (A CT = 1/8, A e = 1). 
The Ising point is easily identifiable since the bound has a "knee" here 4 . After the knee, 
the bound continues to grow linearly and passes via points corresponding to OPEs realized 
in higher minimal models. 

5.3 New stuff 

At this point, one may and should ask: what is the origin of the knee in the D = 2 bound? 
What is so special about A CT = 1/8 that the behavior of the bound changes at this point? 
We would like to show two plots which shed some light on this question. 

The basic idea is to look at the second, after e, scalar in the OPE a x a, which we will 
denote e' . In the free scalar theory, the next scalar in the OPE (5.14) is marginal: 

e' = (dX) 2 , A e , = 2 (free 2D scalar) . (5.15) 

On the other hand, in the 2D Ising model this scalar is strongly irrelevant: 

£' = L_ 2 L_ 2 -1, A £ , =4 (2D Ising). (5.16) 

These facts suggest the following exercise. For each A CT , let us fix A £ to the maximal value 
allowed by the bound in Fig. 3. Let us then ask what is the maximal value of A e / consistent 
with such A CT , A e , and the crossing symmetry of the a four point function. From Eqs. (5.15), 
(5.16) we can suspect that such an upper bound on A £ / should hover around 2 for A CT <; 1/8, 
and then shoot up to about 4. 

The bound on A e > can be easily computed by the algorithm used in [22] to produce the 
bound on A e . Both bounds are plotted in Fig. 4. 5 The bound on A e / completely conforms 

4 In 4D, a similar although less pronounced "knee" has been noticed in a SUSY bound on dimensions of 
real scalars appearing in the chiralxantichiral OPE [28]. It does not correspond to any previously known 
theory and may be a first hint of a new SCFT. 

5 The bounds in Fig. 4 use the linear programming algorithm in the subspace of conformal block deriva- 
tives up to order 10. The best 2D bound on A e presented in [22] used derivatives up to order 12 and was 
somewhat stronger. 
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max A e / (A e — y max) 



max A f 



0.00 



0.05 



0.10 



0.15 



— A 

0.20 CT 



Figure 4: Lower curve: maximal possible value of A e as a function of A CT , computed from the 
crossing symmetry constraint by using the algorithm of [22] . 5 Upper curve: maximal possible 
value of A £ / as a function of A CT and A e (the latter fixed to the maximal value allowed by the 
first bound). The dots are computed; the dashed lines are interpolated. 



with the above expectation: it shows a very steep, essentially step-function, growth around 
A ff = 1/8 from A e , « 2 6 to A e > £ 4. 

A natural interpretation of this plot is that A CT = 1/8 is the dividing line separating 
theories where e is the only relevant operator from theories where necessarily additional 
scalars with A <^ 2 must be present in the a x a OPE. 

Let us test this hypothesis further. We will study the crossing symmetry constraint 
demanding that there should be at most one scalar of dimension < 3 in the a x a OPE. 
In other words, we suppose that A e / > 3. The cutoff value 3 here is picked somewhat 
arbitrarily: it is chosen to exclude the free scalar but to allow the 2D Ising model. 7 As a 
consequence of a very sharp drop in the A e > bound of Fig. 4, the results below will depend 
rather weakly on this value. 

We now ask which region of the (A CT , A e ) plane is consistent with the assumed constraint 
on A e i and the crossing symmetry. Once again, an answer to this question takes only a few 
minutes of your laptop's time to compute via the algorithm of [21, 22]; it is plotted in Fig. 5. 

This plot is interesting in several aspects. In marked difference with Fig. 3, most of the 
allowed region goes away as a result of the A e / constraint. What is left is a curvy triangular 
region localized entirely at A CT ^1/8. To be precise, the tip of the triangle is found located 
at the point A CT ~ 0.124, A £ w 0.996, within 1% from the exact 2D Ising model values (4.1). 

For A CT > 1/8, the upper edge of the allowed region traces the corresponding part of the 
bound in Fig. 3. Moreover, in this range of A CT we also obtain a lower bound on A e . The 
existence of this bound is a consequence of the assumed gap between the dimension of A £ 

6 There is no contradiction between the bound dipping a little below 2 and the fact that A e * = 2 for free 
scalar. In fact A e is fixed to its maximal possible value in this exercise, which is somewhat larger than the 
free scalar line (5.12). 

7 Another banal theory excluded by this constraint is the "generalized free scalar" , i.e. a Gaussian scalar 
field of dimension A CT . In this case the OPE a x a contains operators e = : a 2 : and e' = :(d 2 a)a: of dimensions 
2A CT and 2A CT + 2, respectively. 
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i.o ; 
o.8; 

0.6; 

o.4; 
o.2; 

0.00 0.05 0.10 0.15 0.20 A a 

Figure 5: Shaded: the region of the (A CT , A e ) plane consistent with the assumed constraint 
A e / > 3 and the crossing symmetry in 2D CFT. Computed via the algorithm of [21, 22] with 
derivatives up to order 10. The tip of the allowed region is at the point A CT ~ 0.124, A e ~ 0.996. 




2D Ising model 



and A e /. 8 

To summarize: In the previous section we learned that the knee in the bound of Fig. 3 
was at the 2D Ising model point. This was intriguing but not sufficiently precise to determine 
the a and e dimensions, since the knee was not very sharply defined. In this section we 
instead found that by imposing a requirement that e' be strongly irrelevant, one can get 
a much sharper constraint. The 2D Ising is now at the tip of the allowed region, and the 
dimensions are easier to extract. 

What makes this approach useful is that the results depend only weakly on the assumed 
lower bound on A e /. Thus we can use a very rough estimate for A £ /, like the one coming 
from the first order e-expansion, and still get a good determination of A CT and A £ . 



6 The road to 3D 

Proposal: The 3D Ising model critical exponents can be determined with CFT techniques. 
One should just redo the plots of Fig. 3 and especially Fig. 5 for D = 3. One should then 
look for the knee and for the tip. 

As already mentioned, the lower bound on the e' dimension appropriate to use in the 
3D analogue of Fig. 5 can be inferred from the e-expansion: 

A e / = (4 - e) + e + 0(e 2 ) ^ 4 . 9 (6.1) 

Notice that the operator e' in this case is the S = 4 of Eq. (3.2), which is relevant in the 
UV but becomes irrelevant in the IR. The previous equation refers to its IR dimension. 

On the other hand in the 3D Gaussian scalar theory of the type mentioned in footnote 7, 
the operator e' = :(d 2 a)cr: will have an appreciably smaller dimension 2A CT + 2^3. Thus 

8 Somewhat similarly, Ref. [28] has derived lower bound on the OPE coefficients of protected operators 
in the chiralxchiral OPE as a consequence of a gap between protected and unprotected operators. In their 
case the gap followed from SUSY constraints on the OPE structure. 

9 Review [29] gives an estimate A £ / w 3.84(4) from a variety of theoretical techniques. 
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we should be able to differentiate between this boring theory and the 3D Ising by imposing 
a constraint A 6 / ^3.5. Hopefully, we will get an allowed region in the (A a , A e ) plane which 
will not be very sensitive to this constraint, just as it happened in D = 2. The 3D Ising 
should sit at the tip of this region. 

This project will need some time to carry out, because no compact representations are 
available for the 3D conformal blocks. Such simplified expressions as Eqs. (5.9), (5.11) seem 
to exist only in even dimensions. In D = 3, no particular simplification was found so far, 
and we have to use the following rather more complicated formulas valid for generic D [19] : 10 



g<g(u,v) =u^f«\b,e,A;u,v) 
The scalar function /(°) is given by a double series: 

f°\b,e,A;u,v) 



b = \{A-l\ e = i(A-J). 



(6.2) 



E (A - 6) m (A - e) m (b) m+n (e) m+n _ 
1/ A , -I 1 TV\ I / A "\ U {I V) 



m,n=0 



m\(A + l-\D) m n!(A) 



2m+n 



(6.3) 



where (a) n is the Pochhammer symbol. The higher-spin functions f"\ I > 1 are then 
computed from a two-term recursion relation: 



f\b,e,A;u,v) 
1 A + l-1 



D 
2 



e- 1 



2D-A + 1 



D 



A 



A 



I - 1 



v f {l ~ 1] {b + 1, e + 1, A; u, v) - e + l,A;u,v 



P-i)/ 



e + l-1 
(A + Z- l)(A + /-2 



/^ 1 )(6,e,A; U ^)-^- 1 )(6+l ) e,A;« )W ) 



D 



l)(f -A + e + l-1) (l-l)(D + l-4) 



4 (D - A + I - 2)(D - A + I - 3) 
x uf {l - 2) (b + l,e + l,A;u,v). 



(A-e-l)(e + /-r 



f + l-2)(f + 1-3) 
(6.4) 



Notice that the coefficient of the second term vanishes for I — 1, allowing to start up the 
recursion from just f(°\ 

To apply the algorithm of [21, 22] and compute the bounds, we need to compute Taylor- 
expansion coefficients of the functions around the point u = v = 1/4 (<£>■ z = z = 1/2), 
for any spin and for dimensions allowed by the unitarity constraints (5.6). To do this from 
(6. 3), (6. 4) looks significantly more time-consuming than the corresponding task for the 2D 
and 4D blocks. Nevertheless, it does not look undoable. We hope to turn to this task in 
the immediate future [32]. 

There is however a bright side to Eqs. (6. 3), (6. 4), following from the fact that they are 
valid for any D. The appearance of D as a parameter allows us to analytically continue 
conformal blocks to fractional dimensions. This can be used as a starting point for a 
nonperturbative definition of the Wilson-Fischer fixed points and of field theory in 4 — e 

10 See also [30] for further work on conformal blocks in generic dimensions and in D = 3, and [31] for 
conformal blocks of external operators with spin. 
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dimensions in general. Notice that the crossing symmetry constraint (5.8) is valid in any 
dimension. While our immediate goal is the 3D case, it would be also interesting to analyze 
the D = 4 — e fixed points using our method. As it was noticed in [21], the naive extension of 
the 4D bound to D = 4 — e is in contradiction with the e-expansion; it would be instructive 
to understand how this contradiction gets resolved. 

If our approach is made to work in the simplest case of the 3D Ising model, a general- 
ization to the 3D O(N) model will also be possible. In general, in presence of a continuous 
global symmetry one should classify representations appearing in the OPE. There is more 
freedom since several representations contribute, but there are also more constraints be- 
cause of the symmetry. As a result one can derive a system of equations with the same 
total constraining power as (5.8). Dimensions of each representation (and of the singlet in 
particular) can be bounded individually. The corresponding theory was developed in [26] 
and is valid in any dimension; in [26-28] the idea was shown to work in 4D. 

We argued in this talk that conformal bootstrap methods can give an alternative de- 
termination of the 3D Ising and O(N) critical exponents. While this determination will 
be numerical, its accuracy can be pretty high and improvable with more computer time. 
Importantly, we will be performing a computation which is mathematically well defined, 
since the OPE and a fortiori conformal block decomposition are known to converge [33-35] . 
There will be no need to massage divergent series like in the e-expansion, and there will 
be no related resummation ambiguity. Our short term goal is to try to do better than 
the e-expansion for A CT and A e . The long term goal is to determine the dimensions and 
OPE coefficients of all the operators appearing in the a x a OPE. 11 This would constitute 
essentially a numerical solution of the 3D Ising model at criticality. 
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